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INTRODUCTION 

The design of blades for compressors and turbines ideally requires meth- 
ods for analyzing unsteady, three-dimensional , turbulent viscous flow through 
ui a turbomachine. Clearly, such solutions are impossible at the present time, 

even on the largest and fastest computers. The usual approach at present is 
to analyze only steady flows and to separate inviscid solutions from viscous 
solutions. Three-dimensional inviscid solutions are beginning to be obtained 
with the present generation of computers, but they require excessive computer 
time. So, at present, inviscid analyses usually involve a combination of sev- 
eral two-dimensional solutions on intersecting families of stream surfaces to 
obtain what is called a quasi-three-dimensional solution. 

Since there are several choices of two-dimensional surfaces to analyze and 
many ways of combining them, there are many approaches to obtaining a quasi- 
three-dimensional solution. Most two-dimensional solutions are either on a 
blade-to-blade surface of revolution (Wu's $i surface, ref. 1) or on the 
meridional or midchannel stream surface between two blades (Wu's Sg sur- 
face). However, when three-dimensional effects are most important, significant 
information can often be obtained from a solution on a passage cross sectional 
surface (normal to the flow). This is called a channel solution (fig. 1). 

In the MERIDL program (refs. 4 and 5) a solution to the equations of flow 
on the meridional S 2 surface is carried out. This solution surface is 
chosen when the turbomachine under consideration has significant variation in 
flow properties in the hub-shroud direction, especially when input is needed 
for use in blade-to-blade calculations. The solution can be obtained either 
by the quasi-orthogonal method, which solves the velocity-gradient equation 
from hub to shroud on the meridional stream surface (ref. 2), or by a finite- 
difference method, which solves a finite-difference equation for stream func- 
tion on the same stream surface. The quasi-orthogonal method is efficient in 
many cases and can obtain solutions into the transonic regime. However, there 
is difficulty in obtaining a solution when blade aspect ratios are above 1. 
Difficulties are also encountered with curved passages and low-hub-tip-ratio 
blades. For such cases, the most promising method is the stream-function 
finite-difference solution, which is limited to subsonic flows. 

Stream-function finite-difference programs for flow on the midchannel 
surface of a turbomachine have been reported in the literature. However, many 
of these programs are proprietary or are of limited generality. MERIDL is very 
general and has been thoroughly tested and refined as the result of extensive 
usage at the NASA Lewis Research Center and by industry. 

MERIDL was written for an axial-, mixed-, or radial-flow turbomachine 
blade row, either a compressor or turbine, 0 ^ for an annular duct. Upstream 



and downstream flow conditions can vary from hub to shroud. The solution is 
for compressible, subsonic flow or for incompressible flow. An approximate 
correction for Vjss of stagnation pressure through the blade row is provided. 
The blade row may be either fixed or rotating. The blades may be twisted and 
leaned and can have high aspect ratio and arbitrary thickness distribution. 

The solution obtained by MERIDL also provides the information necessary 
for a more detailed blade-shape analysis o') blade-to-blade surfaces (fig. 1). 
A useful program for this purpose is TSONIC (ref. 3). Information needed to 
prepare all the input for TSONIC is calculated and printed by MERIDL, and is 
stored in the required input format for TSONIC. 


SYMBOLS 

B tangential space between blades, rad 

b stream-channel thickness normal to meridional streamline, meters 
Cn specific heat at constant pressure, J/(kg)(K) 

F vector normal to midchannel stream surface and proportional to tangential 
pressure gradient, N/kg 
m meridional streamline distance, meters 
p pressure, N/meter 2 

R gas constant, J/(kg)(K) 

r radius from axis of rotation, meters 

s distance along orthogonal mesh lines in throughflow direction, meters 
T temperature, K 

t distance along orthogonal mesh lines in direction across flow 
u normalized stream function 

V absolute fluid velocity, meters/sec 

W fluid velocity relative to blade, meters/s'*'' 

w mass flow, kg/sec 

z axial coordinate, meters 

a angle between meridional streamline and axis of rotation, rad; see fig. 4 

e angle between relative velocity vector and meridional plane, rad; see 
fig. 4 

c coefficient in stream-function equation, defined in eq. (3) 
e relative angular coordinate, rad; see fig. 4 
5 coefficient in stream-function equation, defined in eq. (2) 

p density, kg/met er^ 

9 angle between s-distance line and axis of rotation, rad; see fig. 3 
w rotational speed, rad/sec; see fig. 4 

Subscripts : 

m component in direction of meridional streamline 
s component in s-direction 

t component in t-direction 

z component in axial direction 

e component in tangential direction 

Superscripts: 

1 absolute stagnation condition 

" relative stagnation condition 
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METHOD OF ANALYSIS 


Basic Assumptions for MERIDL 

It is desired to determine the flow distribution through a stationary or 
rotating cascade of blades on a midchannel hub-shroud stream surface. The 
following simplifying assumptions are used in deriving the equations and in 
obtaining a solution: 

(1) The flow relative to the blade is steady. 

(2) The fluid is a perfect gas with constant specific heat C p . 

(3) The only forces along a hub-shroud orthogonal mesh line are those due 
to momentum and pressure gradient. 

(4) There is no heat transfer. 

(5) The midchannel surface is a stream surface that has the same shape as 
the blade mean-camber surface, except near the leading and trailing edges, 
where an arbitrary correction is made to match the free-stream flow. 

( 6 ) The velocity varies linearly between blade surfaces. 

(7) The relative stagnation pressure loss is known through the blade row. 

(8) The upstream and downstream boundaries of the solution region are 
orthogonal to the streamlines. 

The flow may be axial, mixed, or radial. Whirl (rV Q ), stagnation pres- 
sure, and stagnation temperature may vary from hub to shroud, both upstream and 
downstream of the blade row. The blade row may be either fixed or rotating, 
with leaned and twisted blades; or there may be no blades at all. Within the 
given assumptions, no terms are omitted from the equations. 

In connection with assumption 3, the viscous forces along a hub-shroud 
orthogonal mesh line are neglected, since these forces are usually very small. 
The viscous forces in the streamwise direction are much larger. The effects 
of these streamwise viscous forces are considered indirectly by specifying a 
streamwise total-pressure-loss distribution. 


MERIDL Stream-Function Solution 


The cylindrical coordinate system is shown in figure 4. However, the so- 
lution is obtained on an orthogonal mesh, as illustrated in figure 3. The 
variables s and t are used to denote distance along the streamwise and 
normal orthogonals, respectively. The stream function satisfies the following 
equation (ref. 4): 
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The derivatives of the stream function satisfy the equations 
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The stream-function equation (1) is a partial differential equation on a 
midchannel hub-shroud stream surface (assumption 5). It is in one unknown (the 
stream function) as a function of two variables and is derived from Wu's mo- 
mentum equation (eq. (96), ref. 1) on what he calls an S 2 surface. Equa- 
tion (1) is nonlinear but can be solved iteratively by the finite-difference 
method when the flow is completely subsonic. 

A finite region (as indicated in fig. ?.) is considered for the solution 
of equation (1). It is assumed that the upstream and downstream boundaries are 
sufficiently far from the blade so as to have a negligible effect on the solu- 
tion. Equation (1) is elliptic for subsonic flow. Therefore, when the flow 
is entirely subsonic, equation (1) can be solved when proper boundary condi- 
tions are specified on the entire boundary of the region. These conditions are 
the values of the stream function or its normal derivative on all four bound- 
aries. The stream function is zero on the hub and 1 at the shroud. On the up- 
stream and aownstream boundaries, the derivative of the stream function normal 
to the boundary is assumed to be zero. This is equivalent to assuming that the 
upstream and downstream boundaries are orthogonal to the streamlines (assump- 
tion 8) . 

The numerical solution of equation (1) is obtained by the finite- 
difference method. A grid must be used for the finite-difference equations. 

The type of grid that was chosen is an orthogonal mesh that is generated by the 
program. The space betwen the hub and the shroud is divided into equal incre- 
ments along several hub-shroud lines. Spline curves are then fit through the 
resulting points to obtain the streamwise orthogonals (fig. 3). The normal 
orthogonals are obtained by a predictor-corrector technique. This technique 
is analogous to the second-order Runge-Kutta method for solving ordinary dif- 
ferential equations, which is known also as the improved Euler method or Heun's 
method (ref. 6). The distances along orthogonal mesh lines are s and t. 

The s-distance is in the streamwise direction and the t-distance is normal to 
this, as indicated in figure 3. With the mesh determined, the finite- 
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difference equations on the orthogonal mesh are given in reference 5. Note 
that the orthogonal mesh is not used as a coordinate system. 

The finite-difference equations are nonlinear since the original equation 
(1) is nonlinear. These equations can be solved iteratively. On the first 
iteration an initial density is assumed; this linearizes some of the terms. 

The remaining nonlinear terms are omitted for the first iteration so that the 
finite-difference equations are entirely linearized. These linearized equa- 
tions are then solved to obtain the first approximate solution forthe stream 
function. This solution provides information that is used to obtain a better 
estimate of the density and an estimate of the other nonlinear terms. The 
equations are then solved again to obtain an improved solution. This process 
is repeated, and by iteration a final converged solution can be obtained if the 
flow is subsonic. 

For each step of this iteration the linearized finite-difference equations 
must be solved. The method used to solve the equations is successive over- 
relaxation (ref. 7) with an optimum overrelaxation factor. Since this is also 
an iterative method, we have two levels of iteration. The overrelaxation is 
performed in the "inner iteration," and the corrections to the nonlinear terms 
are made in the "outer iteration." The inner iteration is internal t« the 
program, so the only iteration apparent to the user is the outer iteration. 

After the stream function is obtained, the velocity distribution is found 
by numerical partial differentiation of the stream function and by using equa- 
tions (5) and (6). The details of the numerical procedure and programming 
technique are described in reference 5. 


Application of MERIDL 

MERIDL can be used both for analysis and as a design tool. When it is 
used for design, other programs should be used with MERIDL. For axial com- 
pressors, reference 8 describes a program that will give blade mean-camber line 
coordinates and thicknesses for an axial compressor blade. This blade design 
can be checked by using the MERIDL program to analyze the flow distribution in 
detail. Usually, changes must be made to the blade design to achieve a desir- 
able flow distribution. These changes may involve more than just the blade 
shape; for example, hub and shroud profile, inlet and outlet whirl (rV«) dis- 
tribution, and loss distribution may have to be changed. Of course, tne accu- 
racy of the MERIDL solution depends on the accuracy of the boundary conditions 
used. 

When a reasonable flow pattern is achieved by the MERIDL flow analysis on 
the midchannel stream surface, more detailed blade surface velocities can be 
obtained by flow analyses on various bl ade-to-bl ade stream surfaces. A useful 
program for this purpose is TSONIC (ref. 3). All of the information required 
to compute input for TSONIC is calculated and printed directly by MERIDL. 
Further changes in blade shape or whirl distribution may be considered at this 
time. Reference 9 (ch. VII) gives information on incidence and deviation for 
good designs and for off-design conditions. 

For cases where the flow is well guided in the channel but has large var- 
iations, both from blade to blade and from hub to. shroud, the CHANEL program 
(ref. 10) is useful. The CHANEL program obtains a solution on a channel cross- 
section surface. The CHANEL program is particularly useful for calculating 
choking mass flow through a blade row. 
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Basic Assumptions for TSONIC 

The simplifying assumptions used are those in references 11 and 12. These, 
assumptions are 

(1) The flow is steady relative to the blade. 

(2) The fluid is a perfect gas (constant C p ) or Is Incompressible. 

(3) The fluid is nonviscous, and there is no heat transfer (therefore, the 
flow is isentropic) . 

(4) The flow is absolutely irrotational . 

(5) The blade-to-blade surface is a surface of revolution. (This does not 
exclude straight infinite cascades.) 

(6) The velocity component normal to the blade-to-blade surface is zero. 

(7) The stagnation temperature is uniform across the inlet. 

(8) The velocity magnitude and direction are uniform across both the up- 
stream and downstream boundaries. 

(9) The only forces are those due to momentum and pressure gradient. 

. The flow may be axial, radial, or mixed and there may be a variation in 
the normal stream-channel thickness b in the through-flow direction. A loss 
in relative stagnation pressure can be specified. 


TSONIC Stream-Function Solution 

The nutation for velocity components is shown in figure 4. For general- 
ity, the meridional streamline distance m is used as an independent coordi- 
nate (see fig. 5). Thus, m and e are the two basic independent variables. 

A stream channel is defined by specifying a meridional streamline radius r 
and a stream-channel thickness b as functions of tn alone (see fig. 6). The 
variables r and b are constant functions of a. 

For the mathematical formulation of the problem the stream function is 
used. The stream function u is normalized so that u is 0 on the upper 

surface of th. lower blade, and 1 on the lower surface of the upper blade. The 

stream function satisfies the following equation (ref. 11). 
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The derivatives of the stream function satisfy 
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If the flow is entirely subsonic, equation (7) is elliptic. Boundary con- 
ditions for the entire boundary ABCDEFGHA of figure 7 will determine a unique 
solution for u. These boundary conditions (ref, 11) are as follows: 
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Boundary segment 

Boundary condition 

AB 

u is 1 less than the value of u on GH 
at the same m coordinate 

BC 

o 

1! 

CD 

u is 1 less than the value of u on EF 
at the same rn coordinate 

DE 
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u is 1 greater than the value of u on 
CD at the same m coordinate 

FG 

u = 1 

GH 
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AB at the same m coordinate 

AH 
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NUMERICAL EXAMPLES 
AXIAL Turbine 

A goal of small-turbine research at the NASA Lewis Research Center is to 
investigate concepts that offer the potential for reducing secondary flows and 
increasing turbine efficiency. One concept that has this potential is stator 
endwall contouring. To provide better understanding of the loss mechanisms 
associated with stator contouring, a program was conducted at Lewis to design 
and evaluate, experimentally and analytically, contoured stator-endwal 1 designs 
for a 1 2 . 77-cm-t ip-diameter axial-flow turbine (ref. 13) as compared with a 
cylindrical stator-endwal! design, The analysis method used MERIDL and TSONIC 
to calculate the blade surface and endwall velocities. Then a boundary layer 
program, BLAVER (ref. 14) was used to calculate stator and rotor displacement 
and momentum thickness, which are then used to calculate profile friction 
losses (including mixing loss) and endwall friction losses. Additional pub- 
lished correlations are used to calculate losses due to incidence, secondary 
flow, rotor tip clearance, disk windage, and exhaust duct friction between 
rotor trailing edge and the downstream measuring station. 

Figure 8 shows a cross section of the turbine, and figure 9 shows the 
stators that were analyzed and tested. The design blade-surface-velocity dis- 
tributions calculated by MERIDL and TSONIC are shown in figure 10. There were 
no static pressure taps on the stator or rotor blades, so only a comparison of 
computed and experimental velocities and flow angles at the stator exit and 
rotor exit could be made. Reference 13 shows very good agreement between ana- 
lytical and experimental velocities and flow angles at the mean section of all 
three stator-rotor configurations. At the hub and tip sections the degree of 
comparison is not as good. This is not unexpected since MERIDL is an inviscid 
program which calculates flow conditions along a hub-to-shroud midchannel 
stream surface. Factors such as endwall boundary layer skewing and three- 
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dimensional tip clearance effects are not accounted for. However, based on the 
overall degree of comparison between the analytical and experimental flow con- 
ditions it was felt that a valid loss assessment could be made. The computed 
efficiencies were slightly higher than their respective experimentally mea- 
sured values. However, the increases in computed stage efficiency for the con- 
toured stators compared to the cylindrical stator was quite close. Although 
there was no comparison of surface velocities or pressures for these stators, 
reference lb' shows extensive experimental comparisons for another stator. 

For practical applications it is important to know the efficiency of a 
turbine over a range of speeds and pressure ratios, not just at the design 
point. The above described analysis has been used to calculate efficiency for 
this turbine at 80 to 100 percent speed and over a range of pressure ratios. 

The predictions are compared with experimental results in figure 11. Further 
details are reported in reference 16. 


Fan Stator 

References 17 and 18 report experimental performance results and analyti- 
cal code calculations for three fan stator designs. These studies were done 
in support of NASA's low-noise, conventional aircraft engine performance. Fig- 
ure 12 shows the flow path for the stage, and figure 13 shows the calculated 
surface velocities for three different stators. Experimental comparisons for 
the surface velocities are not available for these stators. The parameter of 
interest in the calculations was the ratio V maXjSS /V 7 £ as shown in fig- 
ure 13. According to reference 18, this ratio correlated the experimentally 
measured total-loss parameter over a wide range of flows from minimum-loss to 
near-stall operation and for all stators and speeds studied. 
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(b) Contoured stator B. 


Figure 11. - Comparison between measured and predicted stage efficiency for 
turbine stage. 


ORIGINAL PAGE 18 
OF POOR QUALITY 


Instrumentation planes 
(measuring stations); 



Axial distance from rotor-hub leading edge, z, cm 

Figure 12. - Flow path for stages showing axial locations of blading and 'nstru- 
mentation. 
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(a) Stator 9. (b) Stator 9C. 

Figure 13. - Chordwise distribution of surface velocities for S9, 
oints. 
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